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ON THE BL LACERTAE OBJECTS/RADIO QUASARS AND THE FRI/II DICHOTOMY 
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ABSTRACT 

In the frame of unification schemes for radio-loud active galactic nuclei (AGNs), FR I radio galaxies are 
believed to be BL Lacertae (BL Lac) objects with the relativistic jet misaligned to our line of sight, and FR 
II radio galaxies correspond to misaligned radio quasars. The Ledlow-Owen dividing line for FR I/FR II 
dichotomy in the optical absolute magnitude of host galaxy-radio luminosity (M R -L Ra( j) plane can be translated 
to the line in the black hole mass-jet power (Mbh-2jet) plane by using two empirical relations: £2jet-^Rad an d 
Mbh-Mg. We use a sample of radio quasars and BL Lac objects with measured black hole masses to explore 
the relation of the jet power with black hole mass, in which the jet power is estimated from the extended radio 
emission. It is found that the BL Lac objects are clearly separated from radio quasars by the Ledlow & Owen 
FR I/II dividing line in the Mbh-2jet plane. This strongly supports the unification schemes for FR I/BL Lac 
object and FR Il/radio quasar. We find that the Eddington ratios Lboi/^Edd of BL Lac objects are systematically 
lower than those of radio quasars in the sample with a rough division at Lboi/^Edd ~ 0.01, and the distribution 
of Eddington ratios of BL Lac objects/quasars exhibits a bimodal nature, which imply that the accretion mode 
of BL Lac objects may be different from that of radio quasars. 

Subject headings: black hole physics — galaxies: active — galaxies: nuclei — quasars: emission lines — BL Lac- 
ertae objects: general 



1. INTRODUCTION 

FR I radio galaxies (defined by edge-darkened radio struc- 
ture) have lower radio power than FR II galaxies (de- 
fined by edge-brightened radio structure due to compact jet- 
terminating hot spots) (Fanaroff & Riley 1974). Relativis- 
tic jets are observed in many radio-loud active galactic nu- 
clei (AGNs). In the frame of unification schemes of radio- 
loud AGNs, FR I radio galaxies are believed to be misaligned 
BL Lacertae (BL Lac) objects, and FR II radio galaxies cor- 
respond to misaligned radio quasars (see Urry & Padovani 
1995, for a review). Most BL Lac objects have featureless op- 
tical and ultraviolet continuum spectra, and only a small frac- 
tion of BL Lac objects show very weak broad emission lines, 
while quasars usually have strong broad-line emission. The 
broad emission lines of quasars are produced by distant gas 
clouds in broad-line regions (BLR), which are photo-ionized 
by the optical/UV continua radiated from the accretion disks 
surrounding massive black holes. The difference of the broad- 
line emission between radio-loud quasars and BL Lac objects 
may be attributed to their different central engines (e.g., Cav- 
aliere & D'Elia 2002; Cao 2002, 2003). 

The unified scheme of BL Lac objects and FR I radio galax- 
ies have been extensively explored by many previous authors 
with different approaches, such as the comparisons of the 
spectral energy distributions (SEDs) in different wavebands 
(e.g., Owen et al. 1996; Capetti et al. 2000; Bai & Lee 2001), 
the radio morphology, the radio luminosity functions (LFs) 
(e.g., Padovani & Urry 1991; Kollgaard et al. 1992; Laurent- 
Muehleisen et al. 1993) and the optical line emission (e.g., 
Marcha et al. 2005). Padovani & Urry (1992) derived the ra- 
dio LFs of flat spectrum radio quasars (FSRQs) and FR II 
galaxies from a sample of radio-loud AGNs. They consid- 
ered a two-component model, in which the total luminosity 



is the sum of an unbeamed part and a beamed jet luminos- 
ity. The beamed LFs of FR II radio galaxies are consistent 
with the observed LFs of FSRQs and steep spectrum radio 
quasars (SSRQs), which strengthens the unification of FR II 
galaxies and radio quasars (see Padovani & Urry 1992, for 
the details). Similar analyses were carried out on the relation 
between FR I galaxies and BL Lac objects (Padovani & Urry 
1991; Urry & Padovani 1995), which is also consistent with 
the unification of FR Is and BL Lac objects. Even though the 
main observational features of different types of radio-loud 
AGNs can be successfully explained in the frame of the uni- 
fication schemes, some authors have found observations indi- 
cating that the unification may be more complex than usually 
portrayed in these schemes (e.g., Marcha et al. 2005; Landt 
& Bignall 2008). Landt & Bignall (2008) found that a con- 
siderable number of BL Lac objects can be identified with the 
relativistically beamed counterparts of FR II radio galaxies in 
a sample of BL Lac objects selected from the Deep X-ray Ra- 
dio Blazar Survey (DXRBS). 

Ledlow & Owen (1996) found that FR I and FR II ra- 
dio galaxies can be clearly divided in the host galaxy optical 
luminosity-radio luminosity (Mr-Lr^) plane, by a dividing 
line showing that radio power is proportional to the optical lu- 
minosity of the host galaxy. What causes the FR I/FR II divi- 
sion is still unclear, and there are two categories of models to 
explain it: (1) the morphological differences being caused by 
the interaction of the jets with the ambient medium of differ- 
ent physical properties (e.g., Gopal-Krishna & Wiita 2000); 
and/or (2) different intrinsic nuclear properties of accretion 
and jet formation processes (e.g., Baum et al. 1995; Bicknell 
1995; Reynolds et al. 1996; Ghisellini & Celotti 2001; March- 
esini et al. 2004; Hardcastle et al. 2007). Ghisellini & Celotti 
(2001) used the optical luminosity of the host galaxy to esti- 
mate the central black hole masses of FR I/FR II radio galax- 
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ies, and the bolometric luminosity is estimated from the radio 
power of jets in FR I/FR II galaxies. They suggested that most 
FR I radio galaxies are accreting at lower rates compared with 
FR lis, which could correspond to different accretion modes 
in FR I and FR II radio galaxies. If the black hole is spinning 
rapidly, the rotational energy of the black hole is expected 
to be transferred to the jets by the magnetic fields thread- 
ing the holes, namely, the Blandford-Znajek (BZ) mecha- 
nism (Blandford & Znajek 1977). The jet can also be accel- 
erated by the large-scale fields threading the rotating accre- 
tion disk (i.e., the Blandford-Payne (BP) mechanism, Bland- 
ford & Payne 1982). Cao & Rawlings (2004) found that the 
BZ mechanism for rapidly spinning black holes surrounded 
by advection dominated accretion flows (ADAFs) (Narayan 
& Yi 1995) provides insufficient power to explain the jets in 
some 3CR FR I radio galaxies. Wu & Cao (2008) calculated 
the maximal jet power available from ADAFs around Kerr 
black holes as a function of black hole mass with an hybrid 
jet formation model (i.e., BP+BZ mechanism). They found 
that it can roughly reproduce the dividing line of the Ledlow- 
Owen relation for FR I/FR II dichotomy in the black hole 
mass-jet power (Mbh-(2jet) plane with the mass accretion rate 
M ~ 0.01M E dd, if the black hole spin parameter a ~ 0.9-0.99 
is adopted. This accretion rate indicates that FR I and FR 
II galaxies have different accretion modes, supporting the re- 
sults of Ghisellini & Celotti (2001) and suggesting that FR 
I sources are in the ADAF mode. Wu & Cao (2008) 's re- 
sults imply that the black hole spin may play an important 
role in the jet formation at least for FR I radio galaxies (also 
see Sikora et al. 2007, for the discussion of the impact of black 
hole spin on the jet formation in AGNs). 

In this work, we use a sample of BL Lac objects and ra- 
dio quasars with measured radio power, black hole masses, 
and Eddington ratios, to explore the relationship between BL 
Lac objects and radio quasars, and to compare it with the FR 
I/FR II division. The sample and the estimates of black hole 
mass/jet power are described in §2 and §3. We show the re- 
sults in §4, and §5 contains the discussion. The cosmological 
parameters Ho = 70 km s _1 Mpc _1 , Q M = 0.3 and Cl\ = 0.7 have 
been adopted in this work. 

2. THE SAMPLE 

The host galaxies of 132 BL Lac objects have been ob- 
served with the Hubble Space Telescope WFPC2 by Urry 
et al. (2000), among which there are 48 sources with mea- 
sured redshifts and extended radio emission. We add addi- 
tional 18 BL Lac objects compiled in the work of Wu et al. 
(2008) to the Urry et al. (2000)'s sample, which leads to 66 BL 
Lac objects (including 28 low-energy-peaked BL Lac objects 
(LBLs) and 38 high-energy-peaked BL Lac objects (HBLs)) 
with measured redshifts and extended radio emission data for 
our present investigation. We search the literature for the 
emission line data of these sources, and find 44 sources in- 
cluding 23 LBLs and 21 HBLs. We use the luminosity of nar- 
row line [O II] at 3727 to estimate the bolometric luminosity. 
For the sources in which the emission line data of [O II] be- 
ing unavailable, we estimate the [O II] luminosity using other 
narrow emission lines. 

In order to compare the difference between BL Lac ob- 
jects and radio quasars, we need a sample of radio quasars. 
In this work, we adopt the sample of radio quasars compiled 
by Liu et al. (2006), which is selected from the 1 Jy, S4, and 
S5 radio source catalogs. Their sample consists of 146 ra- 
dio quasars including 79 FSRQs (with a2-8GHz < 0.5) and 67 



steep-spectrum radio quasars (SSRQs) (with a2-8GHz > 0.5). 
All quasars in their sample have estimated black hole masses 
and jet power (see Liu et al. 2006, for the details of the quasar 
sample). 

3. THE BLACK HOLE MASS AND JET POWER 

The relation between black hole mass Mbh and host galaxy 
luminosity L K at K-bwA (Eq. 1 in McLure & Dunlop 2004) 
is derived from M^-Mr by using an average color correction 
of R-K = 2.7 for the same cosmology adopted in this paper. 
We convert this relation back to Mbh -Mr as 

log lo (M bh /M ) = -0.50(±0.02)M«-2.75(±0.53), (1) 

to estimate the central black hole masses of BL Lac objects in 
this sample. For a few BL Lac objects, their black hole masses 
can also be estimated from their stellar dispersion velocity a 
with the empirical M^-u relation. It is found that the black 
hole masses of three BL Lac objects estimated with M^-cr 
relation are roughly consistent with those estimated with Eq. 
(1) (see Cao 2004, for the details, and references therein). 

The jet power can be estimated with the relation between jet 
power and radio luminosity proposed by Willott et al. (1999), 

gje, ^ 3 x 10 38 / 3/2 4u5i (W), (2) 

where L extj i5i is the extended radio luminosity at 151 MHz in 
units of 10 28 W Hz" 1 sr" 1 . Willott et al. (1999) have argued 
that the normalization is uncertain and introduced the factor 
/ (1 < / < 20) to account for these uncertainties. This rela- 
tion was proposed for FR II radio galaxies and quasars. Cao 
& Rawlings (2004) compared the power of the jet in M87 (a 
typical FR I radio galaxy) derived with different approaches, 
and found that Eq. (2) may probably be suitable even for FR 
Is (see Cao & Rawlings 2004, for the details, and references 
therein). Following Cao (2003), we adopt this relation to esti- 
mate the power of jets in BL Lac objects, which is believed to 
be a good approximation if BL Lacs can be unified with FR 
Is. 

For most BL Lac objects, their radio/optical continuum 
emission is strongly beamed to us due to their relativistic jets 
and small viewing angles of the jets with respect to the line 
of sight (e.g., Fan & Zhang 2003; Gu et al. 2006). The low- 
frequency radio emission (e.g. 151 MHz) may still be Doppler 
beamed. We therefore use the extended radio emission de- 
tected by VLA to estimate the jet power, as adopted in Cao 
(2003). The observed extended radio emission is Zf-corrected 
to 151 MHz in the rest frame of the source assuming a e = 0.8 
<J V oc (Cassaro et al. 1999). 

We take the black hole masses of radio quasars from Liu 
et al. (2006), which are estimated from the broad line widths 
of H/3, Mg II, or C IV, as well as the line luminosities of these 
lines (see Liu et al. 2006, for the details). In Liu et al. (2006)'s 
work, the jet power is estimated from the extended radio emis- 
sion at 151 MHz with the formula derived by Punsly (2005), 
which is slightly different from Eq. (2) proposed by Willott et 
al. (1999). To be self-consistent, we estimate the jet power of 
quasars in Liu et al. (2006)'s sample from their extended radio 
luminosities with Eq. (2), which is the same as the estimates 
of jet power for BL Lac objects in this work. 

For BL Lac objects, the observed optical continuum emis- 
sion may be dominated by the beamed synchrotron emission 
from the relativistic jets (e.g., Gu et al. 2006). The narrow- 
line regions are believed to be photo-ionized by the radiation 
from the accretion disk, and the narrow-line emission can be 
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used to estimate the bolometric luminosity for BL Lac ob- 
jects. We convert the luminosity of the narrow-line [O II] to 
bolometric luminosity using the relation proposed by Willott 
etal. (1999), 

L bo i = 5 x 10 3 L [OII] W, (3) 

for the BL Lac objects in this sample. For the objects which 
lack [O II] line emission data, we convert the luminosities of 
other narrow lines ([Olll] or H a +[Nll]) to the luminosity of 
[O II] using the ratios suggested by Zirbel & Baum (1995) for 
FR I galaxies. The narrow-line emission data for the BL Lac 
objects are taken from the literature (Sbarufatti et al. 2006; 
Carangelo et al. 2003; Rector et al. 2000; Rector & Stocke 
2001; Stickel et al. 1993; Marcha et al. 1996; Morganti et al. 
1992). We note that Eq. (3) is derived for FR Ils/quasars, 
while ADAFs may be present in these BL Lac objects. The 
SED of an ADAF is significantly different from that of a stan- 
dard thin disk (e.g., Narayan et al. 1995). Nagao et al. (2002) 
calculated the emission of narrow-line regions photo-ionized 
by two different SED templates respectively, i.e., a standard 
thin disk SED template with a bump in UV/soft X-ray bands 
and a hot ADAF SED template described by a power-law con- 
tinuum in hard X-ray bands with an exponential cutoff. They 
found that the narrow-line regions are more efficiently photo- 
ionized by the ADAF SED template than the standard thin 
disk case (see the bottom panel of Fig. 5 in Nagao et al. 2002), 
which implies that the present estimates on the bolometric lu- 
minosity with Eq. (3) may be over-estimated to some extent 
(a factor of ~ 2 - 3 for the narrow-line regions with hydrogen 
column density < 10 20 cm" 2 ). 

For radio quasars, we estimate their bolometric luminosi- 
ties from the total broad-line luminosities Lblr calculated by 
Liu et al. (2006), as the optical continua for most radio-loud 
quasars may probably be contaminated by the beamed emis- 
sion from relativistic jets. Liu et al. (2006) derived a tight 
con-elation: AL A (5100) = 84.3L^ 98 , for the sample of radio- 
quiet AGNs in Kaspi et al. (2000). Given that the luminosity 
of the broad-line Hp corresponds to ~4 per cent of Lblr (see 
Liu et al. 2006, and references therein) and using the relation 
Lbol — 9AL A (5100) (Kaspi et al. 2000), the bolometric lumi- 
nosity can be estimated as Lb i — 30Lblr- 

4. THE RESULTS 

The division between FR I and FR II radio galaxies is 
clearly shown by a line in the plane of total radio luminosity 
and optical luminosity of the host galaxy (Ledlow & Owen 
1996). The optical luminosity of the host galaxy can be con- 
verted to black hole mass Mbh by using the empirical relation 
(1), while the jet power Qjet can be estimated from the radio 
luminosity with relation (2). Thus, the dividing line between 
FR I and II radio galaxies is translated to 

log0jet(ergs s-')= 1.13 logM bb .(M )+33.18+ 1.50 log/, (4) 

in Mbh-2jet plane (see Wu & Cao 2008, for the details), which 
is modified for the cosmology adopted in this paper. In Figure 
1, we plot the relation between the black hole masses Mbh and 
jet power Qjet for radio quasars and BL Lac objects. It is found 
that BL Lac objects can be roughly separated from quasars by 
the FR I/II dividing line. 

The distributions of Eddington ratios for BL Lac objects 
and quasars are plotted in Fig. 2, where only the BL Lac 
objects with measured line emission have been included, be- 
cause the bolometric luminosity is derived from the emission 
lines for these sources. We estimate the statistical significance 



of a possible bimodal distribution of Eddington ratios for BL 
Lac objects and quasars using the KMM algorithm (Ashman 
et al. 1994). The distribution for the entire sample is strongly 
inconsistent with being unimodal (P-value< 0.001), and the 
KMM algorithm separates the entire sample into two groups. 
No significant difference is found in the distributions of black 
hole masses for BL Lac objects and radio quasars. 
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FIG . 1 . — The relation between black hole mass Mbh and jet power Qjct for 
the BL Lac objects and quasars. The open squares and filled squares represent 
FSRQs and SSRQs respectively, while the circles and triangles represent BL 
Lac objects. The filled circles/triangles represent the LBLs/HBLs with mea- 
sured line emission, while the open circles/triangles the LBLs/HBLs without 
measured line emission. The dashed line represents the Ledlow-Owen divid- 
ing line between FR I and FR II radio galaxies given by Eq. (4). 
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FIG. 2. — The distributions of Eddington ratios (Lbol/^Edd) for BL Lac 
objects (dashed line) and quasars (solid line), respectively. 



5. DISCUSSION 
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Figure 1 shows that the FR I/FR II dividing line given 
by Ledlow & Owen (1996) roughly separates the radio-loud 
quasars from BL Lac objects in the Mbh-2j e t plane, which 
strongly supports the FR I/BL Lac objects and FR Il/radio 
quasars unification schemes. This conclusion is independent 
of the value of the uncertainty factor / in Eq. (2). 

We find that only a small fraction of LBLs/quasars are 
above/below the dividing line, which is similar to the FR I/II 
division (see for instance the Fig. 1 in Ledlow & Owen 1996). 
The HBLs have relatively lower jet power than LBLs, and 
only one HBL appears above the dividing line. This means 
that the BL Lacs/quasars and the FR I/II divisions may be 
true only in a statistical sense. The exceptional sources in the 
Mbh-(2jet plane may provide useful clues to investigations on 
the central engines in radio-loud AGNs (e.g., Cao & Rawlings 
2004; Landt & Bignall 2008). 

In Fig. 2, we show that the distributions of Eddington 
ratios for BL Lacs and quasars in our sample exhibit a bi- 
modal nature. The BL Lac objects are roughly seperated from 
the quasars at Lboi/^Edd ~ 0.01, with most BL Lac objects 
having Lboi/^Edd < 0.01 and almost all the quasars having 
iboi/^Edd ^ 0.01. We suggest that this bimodal behavior of 
the distribution may imply different accretion modes in BL 
Lac objects and quasars, and furthermore the transition be- 
tween the accretion states happens at Lboi/^Edd ~ 0.01 ac- 
cording to Fig. 2. Since this is roughly the critical luminos- 



ity above which ADAFs are not possible (e.g., Narayan & Yi 
1995), this suggests that ADAFs are present in BL Lac ob- 
jects and standard thin disks are in quasars. We note that a 
similar explanation is invoked to explain the FR I/II division, 
in which ADAFs would be present in FR I galaxies while stan- 
dard thin disks are in FR II galaxies (e.g., Ghisellini & Celotti 
2001; Wu & Cao 2008). Interestingly enough, Marchesini et 
al. (2004) found a similar bimodal distribution of Eddington 
ratios for a sample of FR I and FR II radio galaxies. 

As discussed in §3, the bolometric luminosities of BL Lac 
objects may be over-estimated, if ADAFs are present in these 
sources. This would strengthen the bimodality in the distribu- 
tion of Eddington ratios of BL Lac objects and quasars. 

The similarity between the division of BL Lac ob- 
jects/quasars and FR I/II found in this Letter strongly sup- 
ports the unification schemes for FR I/BL Lac object and FR 
Il/radio quasar. 
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